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Transient heat transfer in a semitransparent 
radiating layer with boundary convection and 

surface reflections 
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Abstract Surface convection and refractive index effects are examined during transient radiative heating 
or cooling of a grey semitransparent layer with internal absorption, emission and conduction. Each side of 
the layer is exposed to hot or cold radiative surroundings, while each boundary is heated or cooled by 
convection. Emission within the layer and internal reflections depend on the layer refractive index. The 
reflected energy and heat conduction distribute energy across the layer and partially equalize the transient 
temperature distributions. Solutions are given to demonstrate the effect of radiative heating for layers with 
various optical thicknesses, the behavior of a layer heated by radiation on one side and convectively cooled 
on the other, and a layer heated by convection while being cooled by radiation. The numerical method is 
an implicit finite difference procedure with non-uniform space and time increments. The basic method 
developed in earlier work is expanded to include external convection and incident radiation. 
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INTRODUCTION 

Thermal protection coatings and ceramic components 
are being developed for use at high temperatures in 
aircraft and automotive engines. They can be sub- 
jected to transient heating or cooling by a variety of 
external radiation and convection conditions. Some 
of the materials are partially transparent to radiative 
energy. Within the material radiative transport acts in 
combination with heat conduction. Since the materials 
operate at high temperatures there is considerable 
internal emission which is proportional to the layer 
refractive index squared. Ceramic refractive indices 
range from approximately 1.5 to 3 so internal radi- 
ation fluxes can be large. Since these fluxes depend 
strongly on temperature level, accurate instantaneous 
temperature distributions must be calculated during a 
transient numerical solution or the heat flows, and 
hence the solution will become inaccurate as time 
advances. 

For semitransparent materials the transient thermal 
behavior of single and multiple plane layers, cylinders, 
spheres and a square geometry has been examined in 
the literature for a variety of radiative conditions. 
Compared with steady state, transient behavior has 
been studied to a much smaller extent. Two of the 
early transient studies [1,2] are for a semitransparent 
layer bounded by walls at specified temperatures. The 
temperature of one wall is suddenly changed to initiate 
the transient. Index of refraction effects are not 
included so results apply for a gas with a refractive 
index of one. In a more recent work [3] transient 
temperatures in glass were analyzed by using a zonal 
type of method. Computations were made for a variety 


of conditions. Some results are for a layer confined 
between opaque walls with specified temperatures, 
others are for layer assumed to have opaque surfaces 
with temperatures that depend on external conditions. 
A few comments are made on the effects of increasing 
the layer refractive index from 1.5 (typical for glass) 
to 3.0. A cylindrical geometry is analyzed in ref. [4); 
transient cooling characteristics are obtained for the 
cylinder exposed to a cold rarefied environment. A 
multilayered geometry is analyzed in ref. [5], Transient 
cooling of a layer of liquid drops was investigated in 
ref. [6] relative to the design of a radiator to dissipate 
energy in outer space. The layer cools in the cold 
vacuum of space and does not have boundary con- 
vection or incident radiation. A particular solution 
was obtained in ref. [7] for cooling a layer by radiative 
loss in the limit of zero heat conduction. Transient 
radiative cooling of a two-dimensional square 
geometry with a refractive index of one was calculated 
in ref. [8] by a finite difference method. 

The finite difference procedure that is a foundation 
for the present work was developed in refs. [9, 10] . In 
ref. [9] transient cooling was examined for a semi- 
transparent material with a refractive index of one. 
The method was extended in ref. [10] to larger refrac- 
tive indices that provide internal reflections. In refs. 
[9, 10] the layer was initially hot and was cooled by 
exposure to a cold vacuum environment. The present 
work includes radiative and convective heating or 
cooling on each boundary. The layer temperature is 
initially uniform. As the transient begins, the bound- 
ary regions can be heated or cooled rapidly depend- 
ing on the external radiation and/or convection con- 
ditions. To account for large temperature variations 
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NOMENCLATURE 


a 

absorption coefficient of layer [m '] 

.V 

coordinate in direction across layer [m] 

c 

specific heat of radiating medium 
[W s kg 1 K '] 

A 

dimensionless coordinate, v ID. 

D 

thickness of radiating layer [mj 



£, 

h 

. . , E n exponential integral functions 
convective heat transfer coefficient 
[W m : K '] 

convection radiation parameter. 

Greek symbols 

t: M omittance for layer at uniform 

// 

0 

temperature 
time [sj 

k 

n 

N 

<i 

h oi j 

thermal conductivity of layer 

[W m 1 K ! ] 

refractive index of layer 

conduction radiation parameter, 

k/4<j T? D 

heat flux [W m 2 ] 

k ’d 

P 

o 

T 

optical thickness of layer, aD 
density of layer [kg m ’]; surface 
reflectivity 

Stefan Bolt/mann constant 
[W m : K 4 ] 

dimensionless time, (4 aTfpcD)0. 

<■ ! 

dimensionless flux. qja 7? 



</r 

radiative heat flow per unit area and 




time [W m ~] 

Subscripts 

'/.]« 

q l2 external radiative fluxes incident on 

a, b, c. 

d interfaces of layer, Fig. 1 


sides of layer [W m : ] 

g 

gas for convection at boundary 

R 

radiative source term in energy 

i 

initial condition; incoming radiation 


equation 

i 

the Ah A location 

t 

dimensionless temperature, T; T , 

M 

total number of AT grid points 

T 

absolute temperature [KJ 

n 

at the nth time increment 

T\ 

initial temperature of semitransparent 

o 

outgoing radiation 


layer [K] 

s 

surrounding environment 

T m 

integrated mean temperature [K], 

ss 

at steady state 


/n. = TJ T x 

ut 

uniform temperature condition 


temperature of surroundings [K], 

1,2 

external conditions at two sides of 


; s - TJTi 


layer. 


that may develop for some conditions, non-uniform 
increments are used to concentrate grid points near 
the boundaries. Accurate surface temperatures must 
be obtained to prevent errors in the convective heat 
flows and in reradiation from the volume adjacent to 
the boundaries. 

In the solution the radiative effect of the tem- 
perature distribution surrounding each location is 
integrated to obtain its contribution to the local 
absorbed energy within the medium. The transient 
energy equation is solved with an implicit finite differ- 
ence procedure using this local radiative internal 
energy source that depends on the transient tem- 
perature distribution. Internal reflections have a large 
effect on the spatial distribution of absorbed energy. 
The boundary relations developed in earlier work are 
expanded to account for surface reflections, trans- 
mitted external radiation and convection. 

ANALYSIS 

Energy equation for transient moling 

A plane layer with thickness IX Fig. 1, is a gray 
emitting, absorbing and non-scattering material that 
is heat conducting and has a refractive index n 1. 


Initially the layer is at a uniform temperature 7j so 
/(A\0) = 1. It is then placed in surroundings where it 
can receive diffuse radiant energy q Ti and q r2 on each 
side. The layer can also be cooled or heated on each 



Fag. 1. Geometry, boundary conditions and nomenclature 
for plane layer with externally incident radiation and surface 
convection. 
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side by surrounding gas at temperatures 7^ and T &2 
with heat transfer coefficients h x and h 2 . 

The transient energy equation in dimensionless 
form is [1 1] 


where for conditions including incident external radi- 
ation 

A b — (k'd) (4c) 


dt 

dx 


= N 


, o-t 

dX 2 


-R(t) 


(I) 


where R(t) is the gradient of the radiative flux and is 
a function of X and t 


R(t) = n 2 K u t A (X.x) 


A, = 2 p L .£ 3 lK D ) (4d) 


C,(t) = (\-p, d )q Ti +2n 2 p b K D 


i 

l A (X,x)E 2 (K n X)dX 


(4e) 


- -T^ o . h (t)£,(K„A0 + <?., c (T)£ 2 ['c 1 ,(l -AO] 


+ tr k. 


t 4 (X*,x)E x (K„\X*-X\)dX*y (2) 


The q oh and q vx are dimensionless diffuse fluxes that 
are outgoing from the internal sides of the layer 
boundaries (Fig. 1). They each consist of internally 
reflected energy and externally incident radiation 
transmitted through a boundary. 


Boundary conditions 

Boundary conditions are required for radiation and 
for heat conduction coupled with external convection. 
Radiation passes out of the layer from within the 
material; there is no emission at the boundaries which 
are planes without volume, hence there are no radi- 
ation terms in the surface convective boundary con- 
ditions. The conduction-convection boundary con- 
ditions at jt = 0 and v = D for all times are: 


— k 



n 


Mr fI -7X.x = o,0)] 


(3a) 


-k 


dT 
dx v 


= h 2 [T(x = D,0)-T#]. 

n 


(3b) 


The radiation boundary conditions are developed 
in a manner similar to refs. [10, 12]; they are required 
for the time-dependent fluxes </ oh (r) and q ox {x) in 
equation (2). Using the reflectivities on both sides of 
the boundary surfaces, each q 0 is composed of trans- 
mitted and reflected portions: q oh = (1 — p a )<y r , + p b </i.b 
and q ox = ( 1 - p d )< 7 r 2 + p c <7.,c (Fig- 0- The boundaries 
are assumed to be sufficiently rough that all reflections 
are diffuse; values of p are given in ref. [12]. The 
incident fluxes within the layer, q lb and q ix , consist of 
energy leaving the opposite boundary and attenuated 
through the layer, and energy incident at the boundary 
as a result of emission within the layer. These are 
obtained from the radiative flux equation which is the 
integral of equation (2). As detailed in ref. [12], 


Vo.b(T) 


C | (t) 4- A b C 2 (r) 

”T -A b A~ 


(4a) 


c (t) = 


U 2 (t) + ^ c C|(t) 
\-A h A K 


(4b) 


C 2 { T) = (1 -Pd)Qt2 + 2n 2 PcKn 


X 


'I 

t*(X.x)E 2 [K n (\-X)]dX. 

Jo 


(40 


F'or the specific case of black surroundings, q rX — t 4 x 
and q r2 = t 4 2 . 


Numerical solution procedure 

To derive a transient solution method for equation 
(1), the numerical procedure in ref. [10] is used as a 
basis. A brief summary is given here with modi- 
fications for external radiation and convection. As 
shown in ref. [10], trapezoidal integration of dtjdx 
over a small At is used to advance the temperature 
an increment in time to give At = t„ M — /„ % 
(Ax/2)[(dt/dr) n+x +(dt/dz) tl \' The second derivative 
of temperature at t-(-At (at index m+1), is written 
in terms of f(r) (at index n) and At as 
(d 2 tldX 2 )„ M - d 2 AtjdX 2 -yd 2 t n /dX 2 . The radiative 
source at t + At, R( t + At) = /?„+,, is expressed in 
terms of R{x)=R„ by R n . x = R n + (dRjdt)„At. By 
using equation (1) to eliminate dtjdx in the At relation, 
and applying the two preceding relations in this para- 
graph, an equation for At is obtained: 


Ax i 
1 + — 

<dR\ NAx 



2 \ 


(lx 2 



At 


N 


ert 

dx 2 


-R„ 


(5) 


The subscript / specifies the X location (/ = l at A' = 0, 
/ = M at X = 1). Relations are now developed to 
obtain A t(X,) — At, at the X, at t„; the temperatures at 
x „ + 1 are then t n+ x (A',) = r„(A', ) + Ar(A7). Since all terms 
on the right-hand sides of the preceding relations in 
this section are at x corresponding to the index w, this 
subscript is omitted in the following. 

To obtain a solution, relations are needed in equa- 
tion (5) for d 2 jdX 2 at the internal grid points of the 
layer and at the boundaries. For non-uniform 
increment sizes AA', and AX- in the negative and 
positive directions about each A',, the standard second 
derivative discretization is substituted into equation 
(5) to obtain an equation for At, at the interior points 
2 ^ i^ A/-1. 
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NAt 

A/ ( i 4 

AX, (AX; + AX, ) 




4 


NAt 
AX / A A', j 


A/; - 


-U 


2N _ 
(AAV 4 AAV 


AAV 


/VAt 

AAV (AAV +AA', ) 

AAV 4 A X, 

AA'/ A A", 

t, » “ 


A 


4 


AX, 


— R,1 


Ri>- (6) 


To account for convection at each boundary, equation 
(6) is replaced by special forms developed from the 
boundary conditions in equations (3a) and (3b). At 
i= 1, by using the definition of A/, equation (3a) in 
dimensionless form — ctjcX |, t . , = (//,/4/V)(/ g | — /,) is 
written in terms of At by taking the difference in values 
at t4 At and r: 


vAt 
dX , 


- i 


H , 

4 AT 


-A/, 


(7) 


In the finite difference procedure the d 2 AtjdX 2 and 
d 2 t/dX 2 are needed at / = 1 for use in equation (5). 
Expanding for either £ = Ar or r about / = 3/2: 


£c 

3AT 


dc 1 


|J5 

It 

H 7 

ox 2 


- i 


AAV 

i 


4 • ■ - . (8) 


NAt 

(AAV,) 2 


Af A , j 4 


At//)/? 


2 V <9/ 


/VAt 

(AA' V ,) 2 


14 


H_AX x 

AN 


A /, 


— At 


IN 


JAA V)- 


; ( f v, i __ / .v, ) 


// 

4 U a2 t„)-Rs, 

2AX Xf 


i=M. (10b) 


Equations (6) and (10) provide a tridiagonal system 
for obtaining Af h A/ 2 ...,A/ iW ; the coefficients are 
given in the Appendix. 

dRjct is needed for the h, coefficients. From equa- 
tion (2): 



Arrant* (T.r) 


2 | dt 


4; (a.'i>T) 


d^(i) 

dr 


E 2 [ 


+ 2 /iVn 


n , ct(X*. T) 

x /■ (A*,r) 

a’i 

(K D iJf*-Jf|)d^* 

Jo & 

X* 

_ 


Xt(X.x) 


VT 

X _ 


From equation (4) the time derivatives of q ob and <y ot 
in equation (11) are: 


The quantity (^ 2 — Ct)/AA/ is substituted in equation 
(8) for the first derivative at / = 3/2, and either equa- 
tion (7) or equation (3a) is substituted for dAt/dX or 
dtjvX at / = 1. After rearrangement. 


d^o. b (T) _ (dC ,(r)/dT)4v4 b (d C 2 (T)/dT) 

dr 1 -A h A c 

d^ 0 .c(T) _ (dCVlTj/dt) 4 A c (d C \ (T)/dr) 
dr 1 —A b A c 


2 

At z —Ai\ 

H i . 



— A/j 

(.1 _ AA? 

AAV 

AN 1 


(9a) 


? 2 t 

dX 2 


2 

AXt 


h-Jy 

AXt 


4N 


h) 


■ (9b) 


For i = 1 substitute equations (9a) and (9b) into equa- 
tion (5) to give, after rearrangement. 


where 


dC, (t) 
dr 


. r St(X, r) 

t (A,t) — 


E 2 (k u X)AX 


dC,(i) 

dr 




X 


r\ 


J, 


t\X. T) 


SKX.x) 

VT 


E 2 [K D (\-X)]dX. 



H\ AX} 
+ 4.V ” 



NAt 

(Axn 2 


At 2 


= At 


2N 


(AT/ )‘ 


~(h ~ t\ ) + 


H x 


2AX 


t (Vi " ^ i ) — R i 


/ = 1 . (10a) 


Similarly, at / = A/, 


The tridiagonal array is solved using the algorithm in 
refs. [13, 14]. At each X h A r, is added to t, to advance 
the temperature to the next time. 

To evaluate the radiative source term /?(/), equation 
(2), and its derivative dRjdt, equation (11), an accu- 
rate integration method is required. Special treatment 
as described in ref. [6] is used as X * approaches X 
since £ t (0) = f ■ An analytical integration is used for 
a very small region about the singularity and Gaussian 
integration (using an IMSL subroutine) is used start- 
ing at this very small distance away from the singu- 
larity. Values of the functions at the unevenly spaced 
locations in the Gaussian subroutine are found from 
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the grid point values by cubic spline interpolation. By 
trying various AA' values and numbers of grid points 
it was found that 50 increments across the layer gave 
accurate results for the k d considered here. Small 
values of AX = 0.1 were used for 10 increments near 
the boundaries where temperature variations can be 
large. To avoid numerical instabilities that can arise 
for a complex integrodifferential equation as given by 
equations (1) and (2), a small time step At = 0.005 
was used for the first 20 At; after that Ar = 0.01 was 
used. Typical computing times to reach t — 2.0 were 
4 min on a Cray X-MP. In most instances the tem- 
perature distributions for t = 1.5 were less than 1% 
away from steady state. 


Transient energy balance 

The transient temperature distributions were used 
in an overall energy balance to check numerical solu- 
tion accuracy at each time step. The instantaneous 
energy rate incident by radiation and added by 
convection is in dimensionless form, <y r i + < 7 r 2 + 
// i [^ 1 ~/(0,T)] + // : [/ g2 -/(l,T)]. This must equal the 
sum of radiative energy reflected and emitted by the 
layer and the transient rate of energy storage. The 
reflected energy is p a <?n +Pd< 7 r 2 - The instantaneous 
radiative flux leaving through both boundaries is 

U-PhWi.h+O-ArWi* = [(l-Pb)/Pb]tfo.b-0-Pa)</rl] 

+ [( 1 — p c )./pJ[</o c — (1 “ AiWd ■ This was evaluated 
using ij oh and q ox from equation (4). Note that in the 
limit when n = 1 , p h — p c = 0 and A h = A c — 0; the p b 
and p c in the radiative flux calculation can then be 
removed algebraically to avoid singular behavior in 
the numerical calculations. The transient energy stor- 
age rate was obtained from (4/Ai){/ m (T + At) — / m (r)] 
where / m (r) is the instantaneous integrated mean tem- 
perature across the layer. The overall energy balance 
was satisfied within 0.5% throughout the transient 
calculations. 


L imit for infinite thermal conductivity 

Some of the transient temperature distributions that 
will be given are for symmetric heat transfer con- 
ditions on both sides of the layer. For some values 
of the parameters the temperature distributions are 
somewhat uniform during a large portion of the tran- 
sient. In the limit of infinite thermal conductivity the 
distribution is uniform at each time. It is of interest 
to compare this limiting solution with appropriate 
transient results. 

For a gray layer at uniform temperature the layer 
emittance, which is equal to its absorptance, was 
derived in ref. [15] as 


£ U t(«A|,) = a llt («, Kfy) 
= [1 -pM] 


1-2 E,(k d ) 


fbM) 


2 L \ (k'd) 


( 13 ) 


where pjn) is given in ref. [12] from the Fresnel 
relations. For symmetric conditions q ri — q r2 = q rs < 
r gi = r g2 = r g , and /?, = h 2 = h. The transient energy 
balance of radiant absorption and emission, con- 
vection and internal energy storage is 

, D d 'T ut 

«uWr ,-R M <rT„,+h(T f -7 ut ) = pc- ^ . (14) 


The a ut = *: ut and equation (14) is placed in dimen- 
sionless form. The t and / ut variables are separated and 
the result integrated to yield the time corresponding to 
the transient temperature f ul : 


t(C) 



di 


it(</rs ~C 4 ) + 


(15) 


Values for c M (n,K D ) are in refs. [10. 15], Results from 
equation (15) are shown in some of the figures that 
follow and, where applicable, very good agreement is 
obtained with the transient finite difference solutions. 


RESULTS AND DISCUSSION 

The transient temperatures begin from a uniform 
initial temperature. For the results in Figs. 2-5, 7 and 
8, the external radiation and convection conditions are 
symmetric on both sides of the layer so the transient 
temperatures are symmetric and the distributions are 
given for one-half of the layer. Results are given at 
five instances during the transient; the distribution for 
the largest time is at or very close to steady state. 
Some of the results for large time were checked with 
steady state calculations using the computer program 
from ref. [12]; excellent agreement was obtained. 

Figures 2-5 give transient temperatures for a 
layer suddenly subjected on both sides to a higher 
temperature radiative environment so that 
q r] = a T^\ = q r2 = oT* 2 i and to convective cooling. 
This could correspond to a film-cooled ceramic com- 
ponent exposed to a high heat flux from a soot filled 
flame. Figure 2(a)-(c) is for optical thicknesses 
k 0 = 0.5, 2 and 10. For K n = 0.5 the layer is somew hat 
optically thin, and radiation effects are limited by 
small absorption. For K n = 2 radiation can penetrate 
the entire layer, but there is significant absorption; 
this yields close to the maximum radiation effect 
throughout the layer. For k d = 10 the layer is some- 
what optically thick and most absorption of external 
radiation is near the boundaries. During the transient 
the layer is heated on each side by black surroundings 
at a temperature 1.5 times the initial temperature 
(q Tl = q r2 — f.5 4 ). Simultaneously the layer is cooled 
by gas on each side at T g , = 7^ = 0.57’. The 
H ] = H 2 = 1 parameters are such that convection is 
somewhat comparable with radiation. There is mod- 
erate conduction as given by the parameter N = 0.1. 
The solid curves are for a refractive index n = 1, the 
dashed curves are for n = 2. 

In Fig. 2(a) where the layer is rather optically thin, 
external radiation passes readily into the interior to 
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9 I , L i i , i , i ■ I 

0 . 0 0 . 1 0 . 2 0 . 3 0 . 4 0 . 5 

Dimensionless coordinate, X = x/D 


provide a fairly uniform internal heat source, some of 
which is removed by reradiation. Energy is removed 
near the wall by convection interacting with conduc- 
tion. The result is that for small r the interior tem- 
perature remains fairly uniform. As t advances, the 
convection-conduction cooling penetrates the layer 
further and at the steady condition the profile has 
become somewhat parabolic as is typical for a layer 
with a uniform internal heat source. When the refrac- 
tive index is increased from one to two the temperature 
profiles are more uniform. This is a result of increased 
radiative energy transfer across the layer arising from 
internal reflections. It is noted that for small t the 
layer heats more rapidly when n - 2. This is caused 
by an increased layer absorptivity for small k d as 
showm in ref. [15]. 

The optical thickness is increased to x,, = 2 in Fig. 
2(b). It is more difficult for radiation to pass into the 



Fig. 2. Effect of refractive index on transient temperature 
distributions in a layer initially at uniform temperature, after 
exposure to external radiation and convective cooling. 
Parameters: N = 0.1, r/ rl = q r2 = 1.5 4 , //, = //,= 1, 

i - ? p2 — 0.5. n - 1 and 2. (a) Optical thickness = 0.5; 
(b) optical thickness = 2; (c) optica) thickness jc w = 10. 


layer interior, and there is stronger absorption in the 
regions near the boundaries. For small x this produces 
a maximum temperature at about X = 0.1. Heat is 
then conducted into the interior, and to the boundary 
where it is removed by convection. As t increases, the 
interior temperature becomes more uniform as a result 
of conduction with lower radiation penetration; at 
steady state the temperature across most of the 
interior is much more uniform than in Fig. 2(a). 
Increasing n makes the temperature more uniform 
across the layer, but in contrast to Fig. 2(a) the layer 
heats more slowly. The horizontal dot-dash lines are 
for the limit of infinite conductivity from equation 
(15), and are for n — 2. Good agreement is obtained 
w ith the mean temperature levels of the dashed curves 
for n = 2. 

The trends in Fig. 2(b) are accentuated in Fig. 2(c) 
where k*, 5 = 10. For a larger optical density of the 
layer, absorption of incident radiation occurs in 
regions closer to the surface and the temperature 
increases strongly in this region for small r. Energy is 
transported into the interior by conduction and 
internal radiation, and the interior temperature dis- 
tribution becomes uniform as steady state is 
approached. The surface temperature is reduced by 
convective cooling. There is a very significant effect of 
refractive index; the profiles for n = 2 are much more 
uniform. 

Since the profiles have larger variations in shape 
when n ~ 1, results for n — 1 are used in Fig. 3 to 
illustrate further the effect of x- D . Since the N = 0. 1 
parameter contains the layer thickness, /), it is best to 
think of all the results as being for the same D. The 
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Fig. 3. Effect of optical thickness on transient temperature 
distributions during radiant heating and convective cooling 
of a layer with refractive index n = E Parameters: TV = 0.1, 
\] Tl = q r2 = 1.5\ //, = H 2 = 1, / gl = / g2 = 0.5. 


K n is then changed by having materials with different 
absorption coefficients. With this viewpoint all of the 
results are for the same thermal conductivity. For 
x = 0.05 the results demonstrate the effect of increased 
radiant absorption near the boundary as k d increases. 
Near steady state, r « 1.5, increasing the optical den- 
sity makes the temperatures almost uniform except 
near X — 0 within a region of decreasing width as k d 
increases. 

Figure 4(a,b) shows the effect of the conduction 
parameter for k',> = 5 and n — 1 and 2; the previous 
results were for TV = 0. 1 . When TV is increased to one 
in Fig. 4(a) the temperatures are fairly uniform 
throughout the transient. This is accentuated by 
increasing n to two in Fig. 4(b). The results approach 
the limiting horizontal lines from equation (15) for 
TV -► co. For TV — 0.01 and n = 1 there are substantial 
temperature gradients near the boundaries [Fig. 4(a)] 
from the interaction of strong radiant absorption and 
convective cooling. The gradients are substantially 
reduced for n — 2 in Fig. 4(b). 

The effect of surface convection is demonstrated in 
Fig. 5 for — 5. Increasing //, = H 2 to five has a 
very strong effect throughout the transient. Large tem- 
perature gradients exist near the boundaries. Because 
of stronger cooling the steady state temperatures 
(t a 1.5), which represent a balance of radiative heat- 
ing and convective cooling, are significantly reduced. 
Comparing Fig. 5(a) and (b) shows the effect of 
increasing n from one to two; for n = 2 the dis- 
tributions are more uniform. 

For Fig. 6 the conditions at the boundaries have 
been changed. The profiles are not symmetric and are 
shown for X = 0 to 1. The layer is heated on the 



9 I , i . i x i , - i r .J 

0 . 0 0 . 1 0 . 2 0 . 3 0 . 4 0 . 5 

Dimensionless coordinate, X = x/D 



0 . 0 0 . 1 0 . 2 0 . 3 0 . 4 0 , 5 

Dimensionless coordinate, X = x/D 

Fig. 4. Effect of conduction-radiation parameter on the tran- 
sient temperature distributions in a layer during radiative 
heating and convective cooling starting from a uniform 
temperature for x D = 5. Parameters: </ rl = q v2 = 1.5 4 , 

//, = // 2 — I, / gl = = 0.5. (a) Refractive index n = 1; (b) 

refractive index n = 2. 

hot side by a radiative flux equal to that from black 
surroundings at T sl = 1.57V There is no convective 
cooling on the hot side, //, = 0; there is only reradi- 
ation on this side. Cooling occurs on the cold side 
(A" = 1) by convection with H 2 = 1 and T g2 = 0.57V 
and by radiation to black surroundings at T s2 = 0.57V 
These conditions simulate heating a ceramic com- 
ponent by exposure on one side to a soot filled flame, 
and providing convective cooling only on the other 
side. 

In Fig. 6(a) there is a large effect of increased refrac- 
tive index in making the temperature distributions 
somewhat more uniform. At the hot side boundary the 
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Fig. 5. Effect of external convection on transient temperature 
distributions in a layer during radiative heating and con- 
vective cooling starting from a uniform temperature for 
k d = 5. Parameters: /V = 0. L i},\ — </ r: — I 5\ — t g2 = 0.5. 

(a) Refractive index n — I; (b) refractive index n = 2. 


lack of convective cooling provides a zero temperature 
gradient since radiation leaves from within the volume 
and not from the surface itself. In Fig. 6(b) the k d is 
increased to 10 and there is more absorption of inci- 
dent radiation near the hot boundary. Early in the 
transient there is a strong temperature rise near the 
hot boundary, while the temperatures decrease sub- 
stantially at the boundary that is convectively cooled. 
Increasing n is not as effective in equalizing tem- 
peratures as in Fig. 6(a) because radiative transfer 
across the layer as augmented by internal reflections 
is reduced by the increased k d . 

Another type of convection radiation interaction is 



Dimensionless coordinate, X = x/D 



Fig. 6. Effect of refractive index on transient temperature 
distributions in a layer initially at uniform temperature after 
exposure to external radiation on one side without convective 
cooling on that side; cooling is by radiation and convection 
at the other side. Parameters: N — 0.1, = 1.5\ r i2 — 0-5, 

H x = 0. H : — 1, 0.5. n = 1 and 2. (a) Optical thickness 

K'„ = 2; (b) optical thickness K' n = 10. 


illustrated by Figs. 7 and 8 where the layer is heated 
symmetrically by convection of a non-radiating hot 
gas. The radiative environment on both sides is cool, 
such as by having a surrounding film-cooled enclosure. 
At steady state, convective heating is balanced by 
radiative cooling from within the layer. Steady state 
temperatures as a function of are shown in Fig. 
7(a) and (b) for n = 1 and 2. For k d = 0 there is 
no emission by the completely transparent layer, and 
convective heating on both sides with = t^ 2 — 2 
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Dimensionless coordinate, X = x/D 

F'ig. 7. Effect of optical thickness on steady state temperature 
distributions in a layer exposed to large convective heating 
while being cooled by radiation. Parameters: jV — 0. 1, 
? s] = t s2 — 0.2, //, = // 2 - 5, / gl = 2 = 2. (a) Refractive 
index n = 1; (b) refractive index n = 2. 


raises the temperature to / ss = 2 throughout the layer. 
When k d -> oc the layer is opaque and the balance at 
the boundaries of convection, absorption and emis- 
sion provides the uniform equilibrium temperature 
that is shown. For k d between these limits the results 
show the effect of internal radiative cooling on the 
steady state temperature distributions. For each k d 
the maximum temperatures are at the boundaries 
where convective heating is being applied, and the 
minimum temperature from radiative cooling is at the 
centerline. In Fig. 7(a) lor n = 1, the lowest centerline 
temperature is for k d ^ 1 . The lowest surface tem- 
perature is for an opaque layer. For = 20 the tem- 
perature distribution is flat over most of the layer, and 
radiation leaves from regions close to the boundaries. 



Dimensionless coordinate, X = x/D 



Dimensionless coordinate, X = x/D 

Fig. 8. Effect of refractive index on transient temperatures in 
a layer initially at uniform temperature after exposure to 
large external convective heating and radiative cooling. Par- 
ameters: N = 0.1, / sl = / s2 ~ 0.2, //[ = fl 2 = 5, ? gl = 1^2 = 2, 
n = 1 and 2. (a) Optical thickness k d = 0. 1 ; (b) optical thick- 
ness =1; (c) optical thickness k d — 10. ( Continued 
overleaf.) 


The results change somewhat in Fig. 7(b) where n = 2; 
for k d > 0.5 temperatures are more uniform in the 
central portion of the layer. 

For three of the r* d in Fig. 7, transient temperatures 
are shown in Fig. 8 for n ~ l and 2. The application 
of convective heating produces a rapid initial tem- 
perature rise near the boundary. In a few instances 
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Dimensionless coordinate, X = x/D 

Fig. 8 - continued. 


radiative cooling produces temperatures in the central 
portion of the layer that are a little below the initial 
temperature. By comparing Fig. 8(a) and (b) it is 
evident that increasing A.* n and n provides more uni- 
form temperature distributions. The effect of internal 
radiative cooling is substantial, as discussed for steady 
state results in Fig. 7. 


CONCLUSIONS 

Transient solutions are obtained for a semi- 
transparent radiating and conducting layer heated or 
cooled by having each boundary exposed to radiation 
and/or convection. Results for n = 1 and 2 show' the 
effect of the layer refractive index with boundary 
reflections assumed diffuse. A transient implicit finite 
difference formulation developed in previous work 
was extended to incorporate convective and radiative 
boundary conditions. Computing times on a Cray X- 
MP are about 4 min for a complete transient solution. 
The results at steady state were compared with solu- 
tions from a steady state analysis; excellent agreement 
was obtained. 

To illustrate the radiative behavior and types of 
solutions that can be obtained, three types of tran- 
sients are examined. For the first, a layer initially 
at uniform temperature is placed in a hot radiative 
environment while being cooled by convection at both 
boundaries. The transient temperature distributions 
can have rapidly varying shapes when the optical 
thickness is large so that absorption of incident energy 
is concentrated near the boundaries. An important 
effect of refractive index is that internal reflections 
promote the distribution of radiative energy within 
the layer; this makes the transient temperature dis- 
tributions more uniform. For the second type of tran- 
sient, the layer is heated by radiation on one side and 
is cooled by convection and radiation on the other. If 
the optical thickness is less than about 10, internal 
reflections provided by a refractive index of two have 


a substantial effect in equalizing the temperature dis- 
tributions. For the third type of transient a layer in 
cold surroundings is subjected to heating on both sides 
by a hot transparent gas. While being heated the layer 
is losing energy by radiation. The minimum tem- 
perature is at the layer centerline, and the lowest 
steady state temperature is for k d ^ 1 . For a large k d 
the temperatures are uniform over the central portion 
of the layer and the amplitude of the temperature 
variation is decreased. 
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APPENDIX: COEFFICIENTS FOR A TRIDIAGONAL 
SYSTEM 

The tridiagonal system in equations (6) and (10) has the 
following coefficients; 
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